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ABSTRACT 

This  paper  presents  a  summary >  of  the  main  concepts  and  references  of  the  Quantitative  Feedback  Theory 
(QFT).  It  is  a  frequency  domain  engineering  method  to  design  robust  controllers.  It  explicitly  emphasises 
the  use  of  feedback  to  simultaneously  reduce  the  effects  of  model  plant  uncertainty  and  to  satisfy 
performance  specifications.  QFT  highlights  the  trade-off  (quantification)  among  the  simplicity >  of  the 
controller  structure,  the  minimization  of  the  ‘cost  of  feedback’,  the  existing  model  uncertainty  and  the 
achievement  of  the  desired  performance  specifications  at  every  frequency  of  interest.  The  technique  has 
been  successfully  applied  to  control  a  wide  variety  of physical  systems.  After  a  brief  introduction  about  the 
essential  aspects  of  the  QFT  design  methodology’,  including  a  wide  set  of  QFT  references,  this  paper- 
presents  a  new  method  to  extend  the  classical  diagonal  QFT  controller  design  method  for  MIMO  plants 
with  model  uncertainty >  to  a  fully  populated  matrix  controller  design  method.  The  paper  simultaneously 
studies  three  cases:  the  reference  tracking,  the  external  disturbance  rejection  at  plant  input  and  the 
external  disturbance  rejection  at  plant  output.  The  work  ends  showing  several  real-world  examples  where 
the  controllers  have  been  designed  using  QFT  techniques:  an  industrial  SCARA  robot  manipulator,  a 
wastewater  treatment  plant,  a  variable  speed  wind  turbine  of  1.65  MW  and  an  industrial  furnace  of  1  MW. 


1.0  INTRODUCTION 

Much  of  the  current  interest  in  frequency  domain  robust  stability  and  robust  performance  dates  from  the 
original  works  of  H.W.  Bode  (1945)  [1]  and  I.  Horowitz  (1963)  [2],  Since  then,  and  during  the  entire 
second  half  of  the  twentieth  century,  there  has  been  a  tremendous  advance  in  the  state-of-the-art  of  robust 
frequency  domain  methods.  One  of  the  main  techniques,  introduced  by  Prof.  Isaac  Horowitz  in  1959  [24], 
which  characterises  closed  loop  performance  specifications  against  parametric  plant  uncertainty,  mapped 
into  open  loop  design  constraints,  became  known  as  Quantitative  Feedback  Theory  (QFT)  in  the  seventies 
[25-27].  This  paper  presents  a  summary  of  the  main  ideas  and  references  of  the  QFT  methodology. 

The  method  searches  for  a  controller  that  guarantees  the  achievement  of  the  desired  performance 
specifications  for  every  plant  within  the  existing  model  uncertainty.  QFT  highlights  the  trade-off 
( quantification )  among  the  simplicity  of  the  controller  structure,  the  minimization  of  the  ‘cost  of  feedback’ 
(bandwidth),  the  model  uncertainty  (parametric  and  non-parametric )  and  the  achievement  of  the  desired 
performance  specifications  at  every  frequency  of  interest. 

Following  this  introduction,  Section  2  presents  a  brief  description  of  the  essential  aspects  of  the  QFT 
methodology.  Section  3  introduces  a  method  to  design  non-diagonal  QFT  controllers  for  MIMO  systems. 
Afterwards  the  paper  describes  some  real-world  applications  of  the  technique,  carried  out  by  the  author:  an 
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industrial  SCARA  robot  manipulator  in  Section  4,  a  wastewater  treatment  plant  of  5000  mVhour  in 
Section  5,  a  variable  speed  wind  turbine  of  1.65  MW  in  Section  6  and  an  industrial  furnace  of  40  metres 
and  1  MW  in  Section  7.  The  paper  ends  with  a  wide  References  Section  that  includes  a  representative 
collection  of  books  and  papers  related  with  the  theory  and  applications  of  QFT. 

2.0  QUANTITATIVE  FEEDBACK  THEORY 

The  Quantitative  Feedback  Theory  (QFT),  first  introduced  by  Prof.  Isaac  Horowitz  in  1959  [24],  is  an 
engineering  method,  which  explicitly  emphasises  the  use  of  feedback  to  simultaneously  reduce  the  effects 
of  plant  uncertainty  and  satisfy  performance  specifications.  Horowitz’s  work  is  deeply  rooted  in  classical 
frequency  response  analysis  involving  Bode  diagrams,  template  manipulations  and  Nichols  Charts  (NC).  It 
relies  on  the  observation  that  the  feedback  is  needed  principally  when  the  plant  presents  model  uncertainty 
or  when  there  are  uncertain  disturbances  acting  on  the  plant. 

Frequency  domain  specifications  and  desired  time-domain  responses  translated  into  frequency  domain 
tolerances,  lead  to  the  so-called  Horowitz-Sidi  bounds  (or  constraints).  These  bounds  serve  as  a  guide  for 
shaping  the  nominal  loop  transfer  function  L{ s)  =  G(s)  /J(s),  which  involves  the  manipulation  of  gain, 
poles  and  zeros  on  the  controller  G(s).  On  the  whole,  the  QFT  main  objective  is  to  synthesize  (loop-shape) 
a  simple,  low-order  controller  with  minimum  bandwidth,  which  satisfies  the  desired  specifications  and 
tackles  feedback  control  problems  with  robust  performance  objectives. 

In  the  last  few  decades  QFT  has  been  successfully  applied  to  many  control  problems.  A  wide  collection  of 
books  and  papers  about  the  main  aspects  of  the  QFT  methodology,  theory  and  applications,  is  included  in 
the  references  section:  controller  loop-shaping  [41-44],  existence  conditions  for  controllers  [45-47],  multi¬ 
input  multi-output  MIMO  systems  [48-63],  time-delay  systems  [64],  digital  QFT  [65-66],  distributed 
parameter  systems  [67-74],  non-minimum  phase  systems  [75-80],  multi-loop  systems  [81-83],  non-linear 
systems  [84-92],  linear  time  variant  systems  LTV  [93-94],  QFT  software  packages  [95-102],  real-world 
applications  [103,127]. 

A  detailed  study  about  the  history  of  QFT  can  be  found  in  the  papers  written  by  Horowitz  [19-21],  Houpis 
[22]  and  Garcia-Sanz  [23].  In  1992,  Houpis  and  Chandler  organized  in  Wright-Patterson  (Dayton,  Ohio) 
the  first  International  QFT  Symposium  [10].  Since  then,  and  with  the  continuous  support  of  Prof.  Houpis, 
the  Symposia  have  been  organized  every  two  years:  Indiana-USA-1995  [11],  Glasgow-UK-1997  [12], 
Durban  -  South  Africa  -  1999  [13],  Pamplona-Spain-2001  [14],  and  Cape  Town  -  South  Africa  -2003  [15]. 
The  next  one  will  be  in  Kansas-USA-2005. 

To  go  into  the  QFT  theoretical  aspects  in  depth,  check  the  excellent  tutorials  written  by  Horowitz  [16-17] 
and  Houpis  [18].  In  addition,  a  major  analysis  can  be  found  in  the  books  written  by  Horowitz  [3],  Houpis, 
Rasmussen  and  Garcia-Sanz  [4],  Yaniv  [5]  and  Sidi  [6].  Finally,  three  special  issues  of  the  International 
Journal  of  Robust  and  Nonlinear  Control  {Wiley)  describe  some  of  the  more  significant  advances  of  QFT: 
Houpis  1997  [7],  Eitelberg  2001-2002  [8]  and  Garcia-Sanz  2003  [9]. 

The  QFT  design  methodology  is  quite  transparent,  allowing  the  designer  to  see  the  necessary  trade-offs  to 
achieve  the  closed-loop  system  specifications.  The  basic  steps  of  the  procedure  (see  also  Fig.  4)  are 
presented  in  the  following  sub-sections.  They  are: 

•  Plant  model  (with  uncertainty),  Templates  generation  and  nominal  plant  selection .P0(jeo). 

•  Performance  Specifications. 

•  QFT  B  ounds  2?(j  co) . 

•  Loop-shaping  the  controller  G(jeo). 

•  Pre-filter  synthesis  A(jo)). 

•  Simulation  and  Design  Validation. 
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2.1  Plant  Model  and  Templates  Generation 


The  plant  dynamics  to  be  controlled  may  be  described  by  frequency  response  data,  or  by  linear  or 
nonlinear  transfer  functions  with  mixed  (parametric  and  non-parametric)  uncertainty  models.  It  can  be 
defined  taking  into  account  the  parameter  uncertainty  of  the  process  at  every  frequency  of  interest  (coO, 
that  is  to  say  the  plant  uncertainty  templates,  so  that  3P(jcol)={P(jco1),  oi.euQQ. 


The  templates  are  sets  of  complex  numbers  representing  the  frequency  response  of  the  family  of  uncertain 
plants  at  a  fixed  frequency  3/J(jo),),  i.e.  a  template  is  a  projection  of  the  ^-dimensional  parameter  space 
onto  the  Nichols  Chart.  Fig.l  represents  the  QFT-template  of  the  plant  P(s)  =  exp  (-rs)/(s2  +  2^mn  s  +  ) 

with  three  parameters,  two  with  uncertainty  (<^=  0.02,  con  =  [0.7,  1.2],  r=  [0,  2]),  at  eo  =  1  rad/s.  For  more 
information  about  the  QFT  templates  see  [28-34]. 


-350  -300  -250  -200  -150  -100  -50  0 

Phase  (degrees) 


Figure  1 :  Template  of  the  plant  at  co  =  1  rad/s 


2.2  Performance  specifications 

The  standard  two  degree  of  freedom  system  which  best  exemplifies  the  feedback  problem  considered  in 
QFT  is  shown  in  Fig.  2.  It  includes  the  set  of  uncertain  plants,  -3T>(jc0i)={T>(jc0i),  0),euQkj-,  the  loop 
controller  -G-  and  the  pre-filter  -F-,  both  to  be  design,  and  the  sensor  dynamics  -H-.  On  the  other  hand,  R, 
E,  U,  Y  and  N  are  vectors  representing  respectively:  the  reference  input,  the  error  signal,  the  controller 
output,  the  plant  output  and  the  sensor  noise  input.  W,  D\  and  D2  are  the  external  disturbance  inputs.  From 
the  structure  we  can  define  the  Eqs.  (1)  to  (3), 


Figure  2:  Standard  two-degree-of-freedom  feedback  structure 
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Table  I.  Transfer  functions  and  specification  models  (H= =1) 
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To  achieve  reliability  and  robustness,  QFT  deals  with  robust  stability  margins  and  robust  performance 
specifications  (disturbance  rejection,  reference  tracking,  etc)  as  objectives  in  terms  of  the  transfer 
functions  of  Eqs.  (1)  to  (3)  over  the  frequencies  of  interest  (Table  I). 


2.3  QFT-bounds 

For  a  nominal  plant  P0(jco),  member  of  the  family  of  plants  within  the  uncertainty  3P(jc»),  the  QFT 
methodology  converts  closed-loop  system  specifications  and  model  plant  uncertainty  in  a  set  of  constrains 
or  bounds  (Florowitz-Sidi  Bounds)  for  every  frequency  of  interest  that  will  have  to  be  fulfilled  by  the 
nominal  open-loop  transfer  function.  They  are  represented  on  a  Nichols  chart.  Such  a  great  integration  of 
information  in  a  set  of  simple  curves  (the  bounds)  will  allow  designing  the  controller  using  only  a  single 
plant,  the  nominal  plant  P0. 

The  cc>i  plant  template,  3P(jo),)  =  {P(jec>i)},  is  approximated  by  a  finite  set  of  boundary  plants  {Pr(jec>i),  r  = 
l,...,m}.  Each  plant  can  be  expressed  in  its  polar  form  as  Pr(jcOi)  =  p(ec>i)  e1"  (eo;)  =  pZO,  and  likewise  the 
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controller  polar  form  is  G(jec>i)  =  g(©i)  e'r|’  =  gZ<j>.  The  controller  phase  (j)  varies  from  -2 n  to  0.  Therefore, 
for  every  frequency  to,,  the  feedback  specifications  {|7’k(jO)l)|<(5k(o)1),  k=l,...,5}  in  Table  I  -Eqs.  (4)  to  (9)- 
are  translated  into  the  quadratic  inequalities  in  Table  II  -Eqs.  (10)  to  (14)-,  see  [16]. 


k 

Table  II.  Quadratic  inequality  on  g 

Eq. 
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[-41 

l  sfj 
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(12) 

4 

f  2  1  1 

P  --j- 

l 

g2  +  2  p  cos(<|)  +  0)  g  + 1  >  0 

(13) 

5 

f  \  f  \  r  23 

PePd  l~~2  g2+2PePd  Pe  COs(^)  +  Bd )  -  COs(<|>  +  0e )  g+  p2-^f  >0 

l  5sJ  l  §5  J  {  85J 

(14) 

The  format  of  these  quadratic  expressions  is: 

Iai(P’0’Sk’<t)  =  a  g2 +b  g  +  c>  0  (15) 


Chait  and  Yaniv  [35]  developed  an  algorithm  to  compute  the  bounds  based  on  quadratic  inequalities  (see 
Table  II),  simplifying  much  of  the  work  on  traditional  manual  bound  computation.  Taking  these 
inequalities  into  account,  it  is  possible  to  compute  them  at  the  NC.  Once  the  bounds  have  been  calculated 
for  the  performance  specifications,  they  have  to  be  grouped  into  a  single  variable. 

Then,  the  worst  case  bound,  i.e.  the  most  restrictive  one  for  every  phase,  is  computed  for  each  frequency 
of  the  work  array  (see  Figure  3).  For  more  information  about  the  QFT  bounds  see  [35-40]. 

2.4  Controller  design 

In  the  design  stage  (loop-shaping),  the  controller  G( s)  is  synthesized  on  the  NC  by  adding  poles  and  zeros 
until  the  nominal  loop,  defined  as  L0  =  Pq  G,  lies  near  its  bounds.  Loop-shaping  considers  bounds  on  the 
NC  to  express  the  plant  model  with  uncertainty  and  the  performance  specifications  at  every  frequency. 

An  optimal  controller  will  be  obtained  if  it  meets  its  bounds  (over  the  continuous  lines  and  under  the 
dashed  lines  at  every  frequency)  and  it  has  the  minimum  high  frequency  gain  (see  Fig.  3). 

Although  current  CAD  tools  for  QFT  controller  design  are  very  helpful  [95-102],  the  loop-shaping  step 
must  be  still  done  manually  using  designer  skills  and  experience.  Even  keeping  the  controller  structure 
fixed,  automatic  tuning  of  parameters  represents  a  great  challenge  [43-44]. 
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Figure  3:  Loop  shaping 


The  general  formulation  for  the  controller  structure  is  expressed  by  the  following  transfer  function: 


C(s)  =  kG- 


}-CZ 
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vF/l  W  V 


mrp 

s  n 

j= 1 


-+1 


^  mcp  ;  - 

n 

7=1 


(  2  2Re(p/)  ^ 

-5  +  1 


y 


VK7I 


kf 


(16) 


where,  £G  is  the  gain,  Zj  is  a  zero  that  may  be  complex  (ncz,  number  of  complex  zeros)  or  real  (/irz,  number 
of  real  zeros),  and  p j  is  a  pole  (real  or  complex)  with  m,v  the  number  of  real  poles  and  mcp  the  number  of 
complex  poles.  Note  that  the  amount  of  complex  zeros  or  poles  must  be  even,  to  have  pairs  of  complex 
conjugate  numbers  and  obtain  a  polynomial  with  real  coefficients.  Controller  may  have  some  poles  in  the 
origin  and  designer  can  check  the  parameter  r  (usually  0,  1  or  2)  to  set  them.  For  more  information  about 
QFT  loop-shaping  see  [41-44],  and  about  existence  conditions  for  controllers  see  [45-47]. 


2.5  Pre-filter  synthesis 

If  the  feedback  system  involves  tracking  signals,  then  the  best  choice  is  to  use  a  pre-filter  F.  While 
controller  G  reduces  the  uncertainty  and  deals  with  stability,  disturbance  rejection,  etc,  pre-filter  F  is 
designed  to  fulfil  tracking  requirements  [4], 


2.6  Simulation  and  design  validation 

Once  the  controller  design  is  finished,  it  is  necessary  to  analyse  the  behaviour  of  the  system  with  the 
controller  previously  obtained.  Closed-loop  response  at  several  frequencies  and  time  domain  responses 
must  be  checked.  The  analysis  will  be  carried  out  with  the  most  unfavourable  cases  due  to  uncertainty  [4]. 
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Figure  4:  QFT  methodology 
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3.0  QUANTITATIVE  NON-DIAGONAL  COMPENSATOR  DESIGN  FOR  MIMO 
SYSTEMS  [60], [62], 

A  fully  populated  matrix  controller  allows  the  designer  much  more  design  flexibility  to  govern  MIMO 
processes  than  the  classical  diagonal  controller  structure.  This  section  introduces  a  new  methodology  to 
extend  the  classical  diagonal  QFT  controllers  design  for  MIMO  plants  with  model  uncertainty  to  a  fully 
populated  matrix  controller  design.  The  section  simultaneously  studies  three  cases:  the  reference  tracking, 
the  external  disturbance  rejection  at  plant  input  and  the  external  disturbance  rejection  at  plant  output. 
Therefore,  the  role  played  by  the  non-diagonal  controller  elements  gy  (i^j)  is  analysed  in  order  to  state  a 
fully  populated  matrix  controller  design  methodology  for  QFT.  The  definition  of  three  coupling  matrices 
(ciij,  C2ij,  C3y)  and  a  quality  function  <yu  of  the  non-diagonal  elements  come  in  useful  to  quantify  the  amount 
of  loop  interaction  and  to  design  the  non-diagonal  controllers  respectively.  This  yields  a  criterion  to 
propose  a  sequential  design  methodology  of  the  fully  populated  matrix  controller,  in  the  QFT  robust 
control  frame.  As  a  consequence  the  diagonal  elements  gkk  of  the  new  non-diagonal  method  need  less 
bandwidth  than  the  diagonal  elements  of  the  previous  diagonal  methods.  The  work  ends  showing  a  real- 
world  example  (section  4),  where  an  industrial  SCARA  robot  manipulator  is  controlled  using  the  new  non¬ 
diagonal  MIMO  QFT  methodology. 


3.1  Introduction 

Control  of  multivariable  systems  (multiple-input-multiple-output,  MIMO)  with  model  uncertainty  are  still 
one  of  the  hardest  problems  that  the  control  engineer  has  to  face  in  real-world  applications.  Two  of  the 
main  characteristics  that  define  a  MIMO  system  are  the  input  and  output  directionality  -different  vectors  to 
actuate  U  and  to  measure  T-;  and  the  coupling  among  control  loops  -some  outputs  y,  can  be  influenced  by 
several  inputs  Ui,  and  some  inputs  u,  can  influence  several  outputs  y,. 

In  the  last  few  decades  a  very  significant  amount  of  work  in  MIMO  systems,  too  numerous  to  list,  has 
been  done.  Useful  techniques  for  designing  multivariable  feedback  systems  have  been  compiled  in 
excellent  references  written  by  Rosenbrock  [128],  O’Reilly  [129],  Maciejowski  [130],  Skogestad  and 
Postlethwaite  [131],  Houpis,  Rasmussen  and  Garcia-Sanz  [4],  Marlin  [132],  Leithead  and  O’Reilly  [133], 
etc.  Some  of  them  collect  the  frequency  domain  approach,  firstly  introduced  and  adopted  by  Rosenbrock 
and  MacFarlane  in  the  UK  and  by  Horowitz  in  Israel,  in  the  sixties,  when  most  of  the  academic 
community  regarded  the  frequency  domain  as  obsolete  for  MIMO  processes. 

On  the  other  hand,  the  original  works  in  MIMO  systems  were  only  made  for  fixed  plants,  without  any 
uncertainty  in  the  model.  The  first  technique  that  made  a  quantitative  synthesis  and  took  into  account 
quantitative  bounds  on  the  plant  uncertainty,  and  quantitative  tolerances  on  the  acceptable  closed-loop 
system  response,  was  introduced  by  Horowitz  in  the  fifties  [24]  and  subsequently  reinforced  and 
thoughtfully  studied  as  it  has  been  introduced  in  Section  2.  That  technique  is  the  Quantitative  Feedback 
Theory  (QFT). 

Using  MIMO  QFT,  Horowitz  [48-50]  proposed  to  translate  the  original  nxn  MIMO  problem  into  n2 
separate  quantitative  multiple-input-single -output  MISO  problems,  each  with  plant  uncertainty,  external 
disturbances  and  closed-loop  tolerances  derived  from  the  original  problem. 

The  first  improvement  of  the  original  MIMO  QFT  method  was  also  introduced  by  Horowitz  [52].  It 
obtains  a  significant  reduction  of  over  design  in  comparison  with  the  previous  method.  It  considers,  in  the 
successive  steps  of  the  iterative  method,  an  equivalent  plant  that  takes  also  into  account  the  controllers 
designed  in  the  previous  steps.  The  book  by  Houpis,  Rasmussen  and  Garcia-Sanz  [4]  presents  a  detailed 
compilation  of  both  methods. 
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However,  although  such  original  MIMO  QFT  methods  take  the  coupling  among  loops  into  account,  they 
only  propose  the  use  of  a  diagonal  controller  G  to  govern  the  MIMO  plant.  This  structure  can  be  improved 
using  non-diagonal  controllers.  In  fact  a  fully  populated  matrix  controller  allows  the  designer  much  more 
design  flexibility  to  control  MIMO  plants  than  the  classical  diagonal  controller  structure.  The  use  of  the 
non-diagonal  components  can  also  ease  the  diagonal  controller  design  problem.  The  difficulty  is  how  to 
design  the  off-diagonal  elements,  especially  if  one  must  consider  engineering  factors  such  as  cost-benefit 
trade-offs  of  using  cross-feeds,  strong  structure  in  the  plant  uncertainty,  system  integrity  and  plant  input 
signal  levels. 

In  the  last  few  years  some  new  methods  for  non-diagonal  multivariable  QFT  robust  control  system  design 
have  been  introduced.  Again  Horowitz  designed  and  applied  a  procedure  for  non-diagonal  G  controllers 
[51].  The  idea  was  to  insert  a  non-diagonal  matrix  pre-compensator  H  before  the  plant,  to  have  a  new 
effective  plant  Pe  =  P  H.  Later,  Franchek  and  Nwokah  presented  a  sequential  loop  frequency  approach  that 
utilizes  a  fully  populated  matrix  controller  to  meet  performance  specifications,  which  may  include  system 
integrity  requirements  [56,  58].  Boje  utilized  the  Perron-Frobenius  root  interaction  measure  to  design  a 
pre-compensator  that  reduces  the  level  of  coupling  between  loops,  before  a  diagonal  QFT  controller 
matrix  is  attempted  [59].  Yaniv  introduced  an  approach  that  emphasizes  the  bandwidth  of  a  non-diagonal 
pre-controller  multiplied  by  the  classical  diagonal  controller  [57].  Kerr  and  Jayasuriya  presented  a  non¬ 
sequential  MIMO  QFT  methodology  [63]. 

In  this  context,  this  section  goes  on  with  a  previous  work  [60,  62]  and  introduces  a  new  methodology, 
based  on  QFT,  to  extend  the  classical  QFT  diagonal  controller  design  for  MIMO  plants  with  uncertainty  to 
the  fully  populated  matrix  controller  design.  The  work  simultaneously  studies  three  cases:  the  reference 
tracking,  the  external  disturbance  rejection  at  plant  input  and  the  external  disturbance  rejection  at  plant 
output.  It  presents  the  definition  of  three  specific  coupling  matrices  (ciy,  c2,h  c3i j),  one  for  each  case.  They 
come  in  useful  to  quantify  the  amount  of  the  loop  interaction  of  the  system.  Furthermore,  a  quality 
function  qtj  of  the  non-diagonal  elements  gy  (i^j)  for  the  three  problems  is  utilized  to  aid  the  design  of  the 
fully  populated  matrix  controllers.  Based  on  the  above  ideas,  the  work  introduces  a  sequential  design 
methodology  for  non-diagonal  QFT  controllers.  It  contemplates  the  design  of  quantitative  controllers  able 
to  achieve  reference  tracking  and  disturbance  rejection  specifications,  taking  also  into  account  the 
reduction  of  interaction  among  loops.  The  diagonal  elements  gi*  of  the  new  non-diagonal  method  need 
less  bandwidth  than  the  diagonal  elements  of  the  previous  diagonal  methods. 

The  work  begins  with  two  sub-sections  that  formulate  the  coupling  matrices  and  the  coupling  elements  of 
the  control  system.  Then  the  fourth  sub-section  introduces  the  expressions  of  the  optimum  non-diagonal 
controller.  Sub-section  five  analyses  and  compares  the  coupling  effects  among  loops  of  both,  the  classical 
diagonal  method  and  the  new  non-diagonal  methodology.  Sub-section  six  introduces  a  quality  function  to 
quantify  the  amount  of  loop  interaction  and  to  design  the  non-diagonal  controllers.  Sub-section  seven 
presents  the  sequential  procedure  to  design  the  fully  populated  matrix  controller.  Sub-section  eight  deals 
with  some  practical  issues  for  using  the  method.  Section  4  applies  the  new  MIMO  QFT  methodology  to 
control  a  real-world  problem:  a  SCARA  robot  manipulator  [115]. 


3.2  The  Coupling  Matrix 

The  objective  of  this  section  is  to  define  a  measurement  index  (the  coupling  matrix)  that  allows  one  to 
quantify  the  loop  interaction  in  MIMO  control  systems.  Consider  a  nxn  linear  multivariable  system  -see 
Fig.  5-,  composed  of  a  plant  P,  a  fully  populated  matrix  controller  G,  a  pre-filter  F,  a  plant  input 
disturbance  transfer  function  Pdi,  and  a  plant  output  disturbance  transfer  function  Pio,  where  P  e  3P  ,  3P 
is  the  set  of  possible  plants  due  to  uncertainty,  and, 
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The  reference  vector  r’  and  the  external  disturbance  vectors  at  plant  input  d,'  and  plant  output  d„  ’  are  the 
inputs  of  the  system.  The  output  vector  y  is  the  variable  to  be  controlled. 


It  is  denoted 
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as  the  plant  inverse  so  that, 
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(18) 


(19) 


where  A  is  the  diagonal  part  and  B  is  the  balance  of  P  \  and  GA  is  the  diagonal  part  and  6'h  is  the  balance  of 

G. 

The  next  paragraphs  introduce  a  measurement  index  to  quantify  the  loop  interaction  in  the  three  classical 
cases:  reference  tracking,  external  disturbances  at  plant  input,  and  external  disturbances  at  plant  output. 
That  index  is  called  the  coupling  matrix  and,  depending  on  the  case,  shows  three  different  expressions:  C\, 
Ci,  Ci  respectively. 

3.2.1  Tracking 

The  transfer  function  matrix  of  the  controlled  system  for  the  reference  tracking  problem,  without  any 
external  disturbance,  can  be  written  as  shown  in  Eq.  (20), 

y  =  (l  +  P  G)~l  P  G  r  =  Ty/r  r  =  Ty/r  F  r’  (20) 
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Using  Eq.  (18)  and  (19),  Eq.  (20)  can  be  rewritten  as, 

Ty/r  r  =  {l  +  A'1  Gd )~l  A'1  Gd  r  +  (l  +  A'1  Gd )_1  A'1  {Gh  r-{B  +  Gh)  Ty/r  r)  (21) 

In  the  expression  of  the  closed-loop  transfer  function  matrix  of  Eq.  (21),  it  is  possible  to  find  two  different 
terms: 

i.  A  diagonal  term  Ty/r_  d, 

^y/r  d  =  (/  +  A']  Ud)  A'1  Gd  (22) 

that  presents  a  diagonal  structure.  Note  that  it  does  not  depend  on  the  non-diagonal  part  of  the  plant 
inverse  B  nor  on  the  non-diagonal  part  of  the  controller  6'h.  It  is  equivalent  to  n  reference  tracking  SISO 
systems  formed  by  plants  equal  to  the  elements  of  A'1  when  the  n  corresponding  parts  of  a  diagonal  G& 
control  them,  as  shown  in  Fig.  6a. 

ii.  A  non-diagonal  term  Ty/I_b, 

Ty/rb  =  (I  +  A~l  Gd )~l  A'1  [<7b  ~(B  +  Gb)  Ty/r  ]  =  (/  +  A'1  Gd)_1  A'1  C,  (23) 

that  presents  a  non-diagonal  structure.  It  is  equivalent  to  the  same  n  previous  systems  with  internal 
disturbances  q-  r}  at  plant  input  (Fig.  6b). 


r, 


ik 


gii 


Pii 


— ► 


(a) 


n 


Fig.  6  i-th  equivalent  SISO  and  MISO  systems 

In  Eq.  (23),  the  matrix  C\  is  the  only  part  that  depends  on  the  non-diagonal  parts  of  both  the  plant  inverse 
B  and  the  controller  G  b .  Hence,  it  comprises  the  coupling,  and  from  now  on  C\  will  be  the  coupling 
matrix  of  the  equivalent  system  for  reference  tracking  problems, 

C\  =Gb-(B  +  Gb)Ty/r  (24) 
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Each  element  cu j  of  this  matrix  obeys, 

n 

clij  =gij  (1-5ij)-X(p*  +  gik)flg  I1"5*) 
k= 1 

where  5ki  is  the  delta  of  Kronecker  that  is  defined  as, 


J4  =  l<=>k  =  i 
l4i  =0^>k^i 


(25) 


(26) 


3.2.2  Disturbance  rejection  at  plant  input 

The  transfer  matrix  from  the  external  disturbance  at  plant  input  dt  to  the  output  y  can  be  written  as  shown 
in  Eq.  (27), 

}=(l  +  P  Gf'  P  di  =  Tyldi  dt  =  Ty/di  Pdi  d\  (27) 

and  then, 

T,Ul  d,  ={l  +  A-'  Gay  A'  d,-(l  +  A'  Gay  A'  ((B+gJt,^  (28) 

In  that  expression  -Eq.  (28)-  it  is  possible  to  find  two  different  terms: 

i.  A  diagonal  term  Ty/ d;_d, 

r,,dij=(/  +  AJ  Gd('  AJ  (29) 

Again,  Eq.  (29)  is  equivalent  to  n  regulator  MISO  systems,  as  shown  in  Fig.  7a. 

ii.  Non  diagonal  term  Ty/d i  b 

Ty/ii_b  =(l  +  A-lGA  f  A1  (B+  Gb )  Ty/di  =  [l  +  A~l  G& A~l  C2  (30) 

that  presents  a  non-diagonal  structure  which  is  equivalent  to  the  same  n  previous  systems  with  external 
disturbances  c2ij  c/q  at  plant  input,  as  shown  in  Fig.  7b. 

In  Eq.  (30),  the  matrix  C2  comprises  the  coupling,  and  from  now  on  C\  will  be  the  coupling  matrix  of  the 
equivalent  system  for  external  disturbance  rejection  at  plant  input  problems, 

C2=(B  +  Gh)Ty/di  (31) 
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Fig.  7  i-th  equivalent  MISO  systems 


Each  element  C2ij  of  this  matrix  obeys, 

c2ij=E(Ak+gikKj  (i-^ik) 

k=l 

where  5ki  is  the  delta  of  Kronecker  defined  in  Equation  (26). 


(32) 


3.2.3  Disturbance  rejection  at  plant  output 

The  transfer  matrix  from  the  external  disturbance  at  plant  output  d0  to  the  output  y  can  be  written  as 
shown  in  Eq.  (33), 

y  =  (l  +  P  G)  1  da  =Ty/do  da  =  Ty/do  Pdo  da  (33) 

and  then, 

Ty,do  =  (/  +  a->  Gd  )'■  da  +  (/  +  A-'  Gd  )■'  A-'  (b  -  (B  +  Gb )  Ty/do )  d„  (34) 

In  that  expression  -Eq.  (34)-  it  is  possible  to  find  two  different  terms: 
i.  A  diagonal  term  7T/do  d, 

Ty/do  d  =  if  +  ^  Gd  )  (35) 

Once  more,  Eq.  (35)  is  equivalent  to  the  n  regulator  MISO  systems  showed  in  Fig.  8a, 
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ii.  Non  diagonal  term  ry/do_b 

Ty/i0_b  =  (/  +  A'1  Gd )"  A1  [B-{B  +  Gb)  Ty/io ]  =  (/  +  A'1  GA )''  A~l  C3  (36) 

that  presents  a  non-diagonal  structure.  It  is  equivalent  to  the  same  n  previous  systems  with  external 
disturbances  c3y  doi  at  plant  input,  as  shown  Fig.  8b. 


(a) 


n 


Fig.  8  i-th  equivalent  MISO  systems 

In  Eq.  (36),  the  matrix  C3  comprises  the  coupling,  and  from  now  on  it  will  be  the  coupling  matrix  of  the 
equivalent  system  for  external  disturbance  rejection  at  plant  output  problems, 

C,=B-{B  +  Gb)Ty/do  (37) 

Each  element  of  the  coupling  matrix,  c3ij  obeys, 

c3ij  =  Aj  (l-^)-Z  (Pi  +  g; k  )  4j  0  -  3k  )  (38) 

k= 1 

where  5ki  is  the  delta  of  Kronecker  as  defined  in  Equation  (26). 

3.3  The  Coupling  Elements 


In  order  to  design  a  MIMO  controller  with  a  low  coupling  level,  it  is  necessary  to  study  the  influence  of 
every  non-diagonal  element  gy  on  the  coupling  elements  c hj,  c2l}  and  c3y,  defined  in  Eq.  (25),  (32)  and  (38). 

These  elements  can  be  simplified  to  quantify  the  coupling  effects.  Then  it  will  be  possible  to  analyse  the 
loop  decoupling  and  to  state  some  conditions  and  limitations  using  fully  populated  matrix  controllers. 
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To  analyse  the  coupling  elements,  one  Hypothesis  is  stated. 
Hypothesis  HI:  suppose  that  in  Eq.  (25),  (32)  and  (38), 


( *  ^ 

( *  ) 

» 

^ik+gik/hq 

for  k  ^  j,  and  in  the  bandwidth  of  t ; 


(39) 


Note  that  the  above  expression  is  scale  invariant  and  is  typically  fulfilled  once  the  MIMO  system  has  been 
ordered  according  to  appropriate  methods  like  the  Relative  Gain  Analysis  [134],  etc.  Then  the  diagonal 
elements  (y  will  be  much  larger  that  the  non-diagonal  ones  4j, 

1 1  jj  |  »  1 1 kj  | ,  for  k  ^  j,  and  in  the  bandwidth  of  t  y  (40) 

Now,  two  simplifications  are  applied  to  facilitate  the  quantification  of  the  coupling  effects  Ciy,  e2l  j,  C3y. 

Simplification  SI:  Using  the  Hypothesis  HI,  Eqs.  (25),  (32)  and  (38),  which  describe  the  coupling 
elements  in  the  tracking  problem,  disturbance  rejection  at  plant  input  and  disturbance  rejection  at  plant 
output  respectively,  are  rewritten  as  shown  Table  III. 

Simplification  S2:  The  elements  4  are  computed  for  each  case  from  the  equivalent  system  derived  from 
Eqs.  (22),  (29)  and  (35).  The  results  are  shown  in  Table  III. 

Table  III:  Simplifications  to  quantify  the  coupling  effects 


Reference  tracking 

External  disturbances 
at  plant  input 

External  disturbances 
at  plant  output 

Simplification 

SI 

cis=& rh  (pS+^J  : 

(41) 

c2ii=^  (/'* !  g  )  •  i  ui 

(42) 

<>i  P-  -ty  Uj+g,j)  :  i  ••  .i 

(43) 

Simplification 

S2 

l  +  s  jj  Pji 

*-1 

hi=  P>>  *-1  (45> 

l+Si i  P» 

>»=  1  <46> 

,  +  «ii  /hi 

Due  to  Simplifications  SI  and  S2,  the  coupling  effects  c iy,  C2y,  C3y  can  be  computed  as, 
Tracking 


clij  =  gij 


gjj  Ghj  +  gij) 

W  +  Sjj) 


i*j 


(47) 


Disturbance  rejection  at  plant  input 


c2ij 


fej  +  gij) 

(Wj  +  gjj) 


i*j 


(48) 


Disturbance  rejection  at  plant  output 
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c3ij  “  Pi] 


*  (  *  \ 
Pjj  l^ij  +  gijj 

W+2ij) 


i*j 


(49) 


3.4  The  Optimum  Non-diagonal  Controller 


Consider  non-diagonal  controllers  to  reduce  the  coupling  effect,  as  well  as  diagonal  controllers  that  help  to 
achieve  the  loop  performance  specifications.  The  optimum  non-diagonal  controllers  for  the  three  cases 
(tracking  and  disturbance  rejection  at  plant  input  and  output)  can  be  obtained  making  the  loop  interaction 
of  Eqs.  (47),  (48)  and  (49)  equal  to  zero. 

Note  that  both  elements,  p and  /?* ,  of  these  equations  are  uncertain  elements  of  P  .  Every  uncertain 
plant  pj*  can  be  any  plant  represented  by  the  family, 

fe}=  .Pif  (1  +  Aij)  -  0 -|Aij|  -  >  for  i,j  =  l,...,n  (50) 


where  /?,. 


►N 


is  the  nominal  plant,  and  Ap*-  the  maximum  of  the  non-parametric  uncertainty  radii  Aj 


*N  *N 

The  nominal  plants  px j  and  p^ 

follow  the  next  rules: 


that  will  be  chosen  for  the  optimum  non-diagonal  controller  will 


a)  If  the  uncertain  parameters  of  the  plants  show  a  uniform  Probability  Distribution  (Fig.  9a)  -which  is 

*  * 

typical  in  the  QFT  methodology-,  then  the  elements  p^  and  p jj  for  the  optimum  non-diagonal 

*N  *N 

controller  will  be  the  nominal  plants  px j  and  p jj  ,  which  minimise  the  maximum  of  the  non- 

parametric  uncertainty  radii  Ap*-  and  Ap*-  that  comprise  the  plant  templates  (Fig.  9b). 

b)  If  the  uncertain  parameters  of  the  plants  show  a  non-uniform  Probability  Distribution  (Fig.  9c),  then 

*  *  .  .  .  .  *N 

the  elements  px j  and  p  y  for  the  optimum  non-diagonal  controller  will  be  the  nominal  plants  p^ 

*N 

and  p  jj  ,  whose  set  of  parameters  maximize  the  area  of  the  Probability  Distribution  in  the  regions 
[  a-  -  £ ,  aY]  +  8  ]  and  [  <7-  -  s,  a-  +  £  ]  (V  parameter  a^,  A,j,  . . .,  a^,  bjj  . . .)  respectively. 


These  rules  of  selection  will  be  analysed  again  in  Section  4.5,  when  we  compute  the  coupling  effects  with 
the  optimum  non-diagonal  controller. 

Now,  making  Eqs.  (47),  (48)  and  (49)  equal  to  zero  and  using  Eq.  (50),  the  optimum  non-diagonal 
controller  for  each  case  is  obtained. 
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Uniform 

PD 

1  Parameter  a,j 

1  1 

1  1 

1  1 

>  1 

1  1 

1  1 

1  1 

(a) 

3-ij  min  aij  max 

Fig.  9  Probability  Distribution  of  the  parameter  a y,  and  Non-parametric  uncertainty  radii  Ap-}  that  comprise 

the  plant  templates 


3.4.1  Tracking 


gjj 


*N  } 

Pi) 


,N 


Pi) 


,  for  i  *  j 


3.4.2  Disturbance  rejection  at  plant  input 


sT-fp“ 


(  *N^| 

-Pi) 

V  J 


,  for  i  ^  j 


(51) 


(52) 
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3.4.3  Disturbance  rejection  at  plant  output 


opt 

*lj 


=  F 


pd 


Sjj 


*N  A 
Pi J 
*  N 


l  ) 


,  for  i  ^  j 


(53) 


where  the  function  Fpd(A )  means  in  every  case  a  proper  function  made  from  the  dominant  poles  and  zeros 
of  the  expression  A. 


3.5  The  Coupling  Effects 

The  minimum  achievable  coupling  effects  -Eqs.  (54),  (56),  (58)-  can  be  computed  substituting  the 
optimum  controller  of  Eqs.  (51),  (52)  and  (53)  in  the  coupling  expressions  of  Eqs.  (47),  (48)  and  (49) 
respectively,  and  taking  into  account  the  uncertainty  radii  of  Eq.  (50).  Analogously,  the  maximum 
coupling  effect  without  any  non-diagonal  controller  -pure  diagonal  controller  cases-  can  be  computed 
substituting  gjj=0  in  the  Eqs.  (47),  (48)  and  (49)  respectively  -Eqs.  (55),  (57),  (59)-.  That  is  to  say, 


3.5.1 

Tracking 

hi  ».  =  K  (Ajj  -Aij)^jj| 

■&ij  &ij 

(54) 

bL  -o = b  &  -  Ai )- 

(55) 

3.5.2 

Disturbance  rejection  at  plant  input 

C2ij  opt  —  ^ij 

'&1J  &ij 

(56) 

hl,ir,=h(1+A»)| 

(57) 

3.5.3 

Disturbance  rejection  at  plant  output 

bl  =  K  (Ay —  Ajj)  s"jj 

&1J  &ij 

(58) 

b  giro  =  b  (1+ Aij)^jj 

(59) 

where, 


and  the  uncertainty  is:  0  <  |Ay  <  A  p~  ,  0  <  A  J  <  ,  for  i,  j  =  l,...,n 


(60) 
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The  coupling  effects,  calculated  in  the  pure  diagonal  controller  cases,  result  in  three  expressions  (55),  (57) 
and  (59)  that  still  present  a  non-zero  value  when  the  nominal- actual  plant  mismatching  due  to  the 
uncertainty  disappears:  A  -  =  0  and  A  -  =  0  .  However,  the  coupling  effects  obtained  with  the  optimum  non¬ 
diagonal  controllers  -Eqs.  (54),  (56)  and  (58)-  tends  to  zero  when  that  mismatching  disappears. 


3.6  Quality  of  the  Designed  Controller 

Fig.  10  shows  the  appearance  of  three  different  coupling  bands  for  a  common  process.  The  maximum 
|cij|gij=o  -computed  from  Eqs.  (55),  (57)  or  (59)-  and  the  minimum  coupling  effects  without  any  non¬ 
diagonal  controller  gy  limit  the  first  one,  as  well  as  the  second  one  is  bounded  by  the  maximum  and  the 
minimum  coupling  effects  with  a  non-optimum  decoupling  element  gy.  Finally,  the  minimum  coupling 
effect  |  cy  I  gij=gij°pt  with  the  optimum  decoupling  element  g°pt  presents  a  maximum  value,  computed  from 
Eqs.  (54),  (56)  or  (58). 


Fig.  10  Coupling  effect  bands  with  different  non-diagonal  controllers 

From  those  ideas  a  quality  function  qn  is  defined  for  a  non-diagonal  controller  gy  (i^j)  so  that, 


<7y(%)  =  100 


!ogio 

max{ 

h 

gij=<^ 

max{ 

ca 

} 

=  ij  8  ij 

l°gio 

max{ 

cy 

gij=0^ 

max{  c 

\J  ) 

a.^a  .oP*  } 

oij  &ij 

> 


(61) 


The  quality  function  becomes  a  proximity  measure  of  the  coupling  effect  cy  to  the  minimum  achievable 
coupling  effect.  Thus,  the  quality  function  is  useful  to  quantify  the  amount  of  loop  interaction  and  to 
design  the  non-diagonal  controllers.  A  suitable  non-diagonal  controller  will  maximise  the  quality  function 
ofEq.  (61). 


3.7  Design  Methodology 

The  proposed  controller  design  method  is  a  sequential  procedure  closing  loops  [58].  It  uses  a  fully 
populated  matrix  controller  that  does  not  assign  any  special  role  to  the  upper  and  lower  triangular  elements 
of  the  controller  G,  and  in  addition,  it  can  be  used  to  design  the  feedback  controller  of  a  2DOF  structure 
(see  Fig.  5).  First  it  is  necessary  to  fulfil  the  Hypothesis  HI  and  two  new  Hypotheses  H2  and  H3. 
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Hypothesis  H2:  The  plant  P  and  its  inverse  P*  have  to  be  stable  and  do  not  have  any  hidden  unstable 
mode.  Despite  looking  very  restrictive,  note  that  it  is  only  a  sufficient  condition  to  guarantee  the  stability 
of  the  system.  Consequently,  the  designer  must  pay  close  attention  to  systems  with  non-minimum  phase  or 
unstable  elements  [58,  61]. 

In  the  last  few  years,  several  works  have  studied  the  stability  problem  of  MIMO  systems  with  uncertainty, 
using  inverse  plants  in  the  QFT  methodology.  A  deep  analysis  of  the  subject  can  be  found  in  the  next  four 
references:  [61,  106,  4,  5].  The  second  paper  proofs  that  it  is  necessary  and  sufficient  that  the  plant  of  each 
successive  loop  is  stabilised.  The  third  and  fourth  references  expand  the  analysis. 


Hypothesis  H3:  The  plant  P  is  not  'ill-conditioned'  for  any  of  the  possible  plants  in  the  whole  set  3P.  This 
will  guarantee  the  robustness  of  the  design.  It  is  known  [135]  that  'ill-conditioned'  plants,  with  large 
elements  of  the  RGA  -Relative  Gain  Analysis  [134]-  matrix  are  difficult  to  control. 


Methodology’ 


Step  A.  First,  the  methodology  begins  paring  inputs  and  outputs  with  the  RGA  technique  [134],  and 

(  -j-  \  ]  /  j|{  \  ] 

arranging  the  matrix  P  so  that  (y?nJ  has  the  smallest  phase  margin  frequency,  (p22  )  the  next  smallest 
phase  margin  frequency,  and  so  on  [4],  Later,  the  sequential  controller  design  technique  (Fig.  11), 
composed  of  n  stages  -n  loops-,  will  follow  the  next  steps  (B  and  C)  for  every  column  k  =  1  to  n. 


L§nl  0 


...  0 

...  0“ 

§11  §12>  -  0 

...  0“ 

§11 

§12 

§  lk 

...  0 

...  0 

§21  (§22)jj  ■"  ® 

...  0 

§21 

§22 

...  g2k  ... 

=> 

/  ... 

=>  ...  => 

...  0 

...  0 

§kl  §k2/  -  0 

...  0 

§kl 

§k2 

§kk 

...  0 

...  0 

_§nl  §  n2  -  0 

...  0_ 

_§  nl 

§  n2 

§  nk 

and  C  ia 

to  Cln 

Step  B2  and  C2a 

to  C2n 

Step  Bn 

and  Cna  to  C 

r 


&ln 

§2n 


’kn 


Fig.  11  Steps  for  controllers  design 


Step  B.  Design  the  diagonal  controller  gkk.  This  design  stage  of  gkk  is  calculated  through  standard  QFT 

loop-shaping  for  the  inverse  of  the  equivalent  plant  (gkk )  in  order  to  achieve  robust  stability  and  robust 
performance  specifications.  The  equivalent  plant  satisfies  the  next  recursive  relationship  [58], 


bn  Jk  =bnLi - — — - — - - - ; 

L^(i-l)(i-l)  Jk_!  +  L  S"(i-l)(i-l)  J  *_/ 


(62) 


which  is  an  extension  for  the  non-diagonal  case  of  the  recursive  expression  proposed  by  Horowitz  [52]  as 
the  Improved  design  technique,  also  called  Second  method  by  Houpis  et  al.  [4]. 


If  the  control  system  requires  tracking  specifications  as  ^((q)  <  t^1  (jco)  <  &;;(&>)  then,  because 


tY  -^rii  +  Qlii  -Eq.(21)-,  the  tracking  bounds  bn  and  au  will  have  to  be  corrected  with  the  coupling 
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specification  rc m,  so  that: 

b[[  b\\  -  Zclii  ,  Cl[[  “t"  Zclii 

'clii  =  Wii  Clii  ^  rclii 

These  are  the  same  corrections  proposed  by  Horowitz  [48,  52]  and  Houpis,  Rasmussen  and  Garcia-Sanz 
[4]  for  the  original  MIMO  QFT  methods  1  and  2. 

However,  with  the  new  non-diagonal  method  these  corrections  will  be  less  demanding.  The  coupling 
expression  tclii  =  ciu  is  now  minor  than  in  the  previous  diagonal  methods  -compare  Eqs.  (54)  and  (55)-. 
The  off-diagonal  elements  gy  (i^j)  of  the  matrix  controller  will  attenuate  or  cancel  that  cross  coupling. 
Then  the  diagonal  elements  of  the  non-diagonal  method  will  need  less  bandwidth  than  the  diagonal 
elements  of  the  previous  diagonal  methods. 

Step  C.  Design  the  (n- 1)  non-diagonal  elements  glk  (i  ^  k,  i  =  1,2, ...n)  of  the  k- th  controller  column, 
minimising  the  coupling  cik  -computed  from  Eqs.  (47),  (48)  and  (49)-.  To  achieve  this  goal,  the  nominal 
optimum  controller  -Eqs.  (51),  (52)  and  (53)-  must  be  taken  into  account. 

Step  D.  Finally,  the  design  of  the  pre-filter  F  does  not  present  any  difficulty,  if  the  complementary 
sensitivity  function  shows  a  low  level  of  loop  interaction.  Therefore,  the  pre-filter  F  can  be  diagonal. 


(63) 

(64) 

(65) 


3.8  Some  Practical  Issues 

To  use  the  above  controllers  there  is  a  relevant  condition  to  take  into  account:  the  plant  P  and  its  inverse 
P  have  to  be  stable  and  do  not  have  any  hidden  unstable  mode  -The  H2  Hypothesis-  [58].  This  is  only  a 
sufficient  condition.  Hence,  in  some  cases  of  not  fulfilling  this  condition  it  is  even  possible  to  use  a  non¬ 
diagonal  controller  element  like  the  above-mentioned  [63].  In  these  situations  some  difficulties  might 
appear  during  the  sequential  procedure  such  as  the  introduction  of  additional  non-minimum  phase  (nmp) 
zeros  due  to  hidden  unstable  modes  of  the  plant  or  due  to  loop  closures.  This  typically  occurs  in  highly 
coupled  systems. 

To  avoid  this  problem,  the  elements  of  the  resulting  plants  (py[)  must  be  checked  in  every  step  of  the 
design  methodology  to  ensure  that  right  half  plane  transmission  zeros  or  unstable  modes  have  not  been 
introduced  by  the  new  controller  elements  gkk  or  gik,  which  would  obviously  cause  an  unnecessary  loss  of 
control  performance  [58].  If  these  nmp  zeros  appear  due  to  the  designed  controller  elements, 
supplementary  constraints  in  the  determinant  of  PG  should  be  imposed  to  re-calculate  the  controller.  On 
the  other  hand,  if  the  plant  elements,  pi±  or  pik,  are  the  cause  of  the  introduction  of  non-minimum  phase 

elements  in  the  equivalent  plant  (gk|)  ,  the  theory  proposed  for  nmp  MISO  feedback  systems  by 
Horowitz  et  al.  [3,  75-79]  and  modified  by  Chen  and  Balance  [80]  can  be  applied  to  properly  design  the 
controllers  in  the  loop  shaping  step. 

At  the  same  time,  arbitrarily  picking  the  wrong  order  of  the  loops  to  be  designed  can  result  in  the  non¬ 
existence  of  a  solution.  This  may  occur  if  the  solution  process  is  based  on  satisfying  an  upper  limit  of  the 
phase  margin  frequency  eo<|),  for  each  loop.  Hence,  Loop  i  having  the  smallest  phase  margin  frequency  will 
have  to  be  chosen  as  the  first  loop  to  be  designed.  The  loop  that  has  the  next  smallest  phase  margin 
frequency  will  be  next,  and  so  on  [4], 
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It  is  important  to  notice  that  the  calculation  of  the  equivalent  plant  )  usually  introduces  some  exact 
pole-zero  cancellations.  That  operation  can  be  precisely  done  using  symbolic  mathematic  tools,  but  could 
be  erroneously  done  when  using  numerical  calculus  due  to  the  typical  round  errors. 

Finally,  every  step  of  the  proposed  methodology  can  be  aided  by  existing  QFT  software  packages  [95- 

102]. 

4.0  NON-DIAGONAL  MIMO  QFT  CONTROLLER  FOR  A  SCARA  ROBOT  [62] 


The  non-diagonal  MIMO  QFT  controller  design  technique  is  applied  to  control  a  real-world  problem:  a 
SCARA  robot  [115],  Fig.  12  shows  the  AdeptOne  robot  manipulator,  and  the  two  joint  angles  -Si  and  82- 
that  are  considered  in  this  example.  In  order  to  present  the  control  of  two  joints  of  the  SCARA  robot,  the 
plant  model  and  the  desired  performance  specifications  are  introduced  in  the  next  sections. 


Fig.  12  Adept  One  SCARA  Robot 


4.1  Plant  Model 


The  Lagrange  equations'  method  is  used  to  find  Eqs.  (66)  and  (67),  which  describe  the  non-linear  dynamic 
behaviour  of  the  two-link  system  [115],  The  real  inputs  are  the  torques  i]  and  X2  -applied  through  power 
amplifier  as  Ui  and  u2-  commanded  by  electrical  motors  on  joints  1  and  2,  and  the  outputs  are  angles  81  and 
82. 


[aL  +2  a3  cos(52)]  8j  +[a2  +a3  cos(52)]  S2  +  Vj  8j  +1^  sgr^Sj)  =  t1  =  —  (66) 

k 
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[a7  +  a3  cos (8,)]  Sj  +a2  57  +  v7  SL  +  p2  sgn(82)  =  x2  =  — 

' . '  "  “  .  k 


(67) 


where  k  is  the  power  amplifiers'  gain,  v;  are  coefficients  of  viscous  friction,  p,  Coulomb  friction 
parameters  associated  with  link  i,  and, 


al  =I1  +I2  +ml  x\  +m2  W  +xl) 
2 


a2  =I2  +m2  x2 

a3  =  m2  ^1  x2 


(68) 


denoting  I;,  m,  and  x,  as  the  moment  of  inertia,  mass  and  position  of  the  i-th  link  respectively,  and  f  as  the 
length  of  link  1 . 

Input  signals  Ui  and  u2  will  be  computed  in  counts  [ct]  and  will  be  commanded  to  the  robot  motors  by  the 
amplifiers.  After  a  robust  parameter  identification  process  [115],  the  coefficients  of  the  robot  model,  with 
a  uniform  Probability  Distribution,  were  found  (Table  IV). 

Table  IV:  Coefficients  of  uncertain  plant.  Uniform  Probability  Distribution. 


Minimum 

Maximum 

Nominal 

ark  [ct-s2/rd] 

719 

813 

766 

a2-k  [ct-s2/rd] 

186 

200 

193 

a3-k  [ct-s2/rd] 

134 

230 

182 

Vrk  [ct-s/rd] 

67 

381 

224 

v2-k  [ct-s/rd] 

11.6 

91.9 

51.75 

Pi-k  [ct] 

344 

358 

351 

ft2-k  [ct] 

262 

323 

292.5 

Now  it  is  possible  to  consider  the  Coulomb  frictions  as  disturbances  and  the  cosine  value  of  S2  as  an 
uncertain  parameter  h  between  -1  and  +1.  Taking  into  account  Eqs.  (66)  and  (67),  it  is  easy  to  find  the 
following  linear  transfer  functions,  which  are  the  elements  of  the  plant  P  defined  as, 


"*i' 

=  P 

ui 

_ 

Pn 

Pn 

U1 

_U2  _ 

_Pn 

P  22  _ 

_U2  _ 

,  \  a7  s  +  v7  1 

Pii(s)  =  TaG0"  k 


P\2  00  = 

/J2l(s)  =  - 

^22  (s)  = 


'  (a2  +  a3  h)  I 

A(s)  k 

-(a2  +a3  h)  | 

A(s)  k 

(a!  +  2  a3  h)s  +  V[  1 
s  A(s)  k 


(69) 

(70) 

(71) 

(72) 

(73) 
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where, 


A  =  t]2  s2+m  s  +  rio 

with  the  following  coefficients, 


(74) 


r|2  =  a2  (ai  +2  a3  h)-(a2  +a3  h)2 
rii  =a2  Vi  +v2  (aj  +2  a3  h) 
no  =vj  v2 


(75) 


4.2  Performance  Specifications 


The  desired  performance  specifications  for  the  SCARA  robot  manipulator  are  the  following, 

i.  Robust  Stability.  |/i;(j®)|  <  1.2  ,  for  i  =  1,2  V®  ,  involves  a  lower  phase  margin  of  at  least  50°  and  a 
lower  gain  margin  of  at  least  1.833  (5.26  dB). 

ii.  Control  effort  constraint.  Control  signals  have  to  be  lower  than  32767  [ct]  for  disturbance  rejection  at 
plant  output  of  about  20°. 

Hi.  Disturbance  rejection  at  plant  input.  The  maximum  allowed  error  has  to  be  of  30°  for  torque 
disturbances  of  1000  [ct], 

iv.  Loop  Coupling.  Reduction  of  coupling  effect  as  much  as  possible. 

v.  Tracking  specifications.  |7’y/r(j®)|  has  to  achieve  tracking  tolerances  defined  by, 


«ii  (®)  ^  kf  0  ®)|  ^  h  (®)  for  i  =  1 , 2 


(76) 


where, 


12.25  [(j®)/30  +  l] 
(j®)2  +5.25  (j ®)  + 12.25 


«ii  (®)  = 


2.25 


[(j®)2  +4.5  ( j ®)  +  2.25]  [( j ®)/10  + 1] 


(77) 

(78) 


The  above  specifications  are  limited  by  the  achieved  sampling  time  for  the  practical  implementation, 
which  is  actually  10  ms. 


4.3  Controller  Design 


•  Step  A:  Coupling  analysis  and  pairing 

The  first  step  is  the  RGA.  This  analysis  yields  a  very  obvious  result:  angle  81  will  be  controlled  by  motor  1 
(iq),  and  angle  S2  by  motor  2  ( u2 ).  The  first  element  A,n  plotted  in  Fig.  13  also  shows  that  the  robot  arm 
presents  a  very  coupled  behaviour.  At  low  frequencies  -below  0.06  rad/s-  the  coupled  behaviour  is  very 
low,  but  as  far  as  the  frequency  increases  the  system  presents  a  more  coupled  dynamics.  The  required 
bandwidth  of  the  system  derived  from  tracking  specifications  lies  between  approximately  2  and  3.5  rad/s. 
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In  those  frequencies  the  maximum  value  of  Xu  is  greater  than  4.5. 


Frequency  (rd/s) 


Fig.  13  Element  Xu  of  the  Relative  Gain  Analysis  Matrix 

•  Step  B.l:  Design  of  the  first  loop  controller,  gn(s). 


Taking  into  account  the  model  uncertainty  of  (/?*, )  'and  the  desired  specifications  for  the  Loop  1,  the 
controller  of  Eq.  (61)  is  found,  satisfying  all  the  performance  specifications  -see  Fig.  14a-. 


fTnOO 


1,65  s2  +4.3840  s  +  2.6190  9 

(s2  +829.2  s  + 1.545  105)s 


(79) 


•  Step  C.l:  Design  of  the  decoupling  element  of  control  effort  ui  on  angle  82. 

Taking  into  account  the  optimum  controller  for  reference  tracking  problems  of  Eq.  (51),  the  controller  g2i 
of  Eq.  (80)  is  designed  minimising  the  coupling  effect  c2i.  Fig  14c  shows  the  frequency  plot  of  the 
obtained  coupling  reduction. 


.  .  3860  s2  +10300  s  +  6130 

g2i(s)  = - 2 — ; - 5 - 

848  10~6  s3  +0.00703  s2 +1.31  s  +  0.346 
•  Step  B.2:  Design  of  the  second  loop  controller,  g22(s). 

The  following  equivalent  plant  (p22 )  *  derived  from  Eq.  (62)  is  calculated, 


r  *e  1  r  *  1  ([p2i]i+k2i]i)([pi2]i+ki2]i) 

Lp22J2  -  Lt22Ji  r— j 

[puJi+bnli 


(80) 


(81) 


The  controller  of  Eq.  (82)  is  found,  satisfying  all  the  performance  specifications  -see  Fig.  14d-. 

88.1s2  +225  s  +  110 

S22  (s)  =  /  ,  , - ~ - v-  (82) 

(l0“6  s2 +0.371  10~3  s  + 0.0344) s 

•  Step  C.2:  Design  of  the  decoupling  element  of  control  effort  m2  on  angle  8]. 

The  controller  gn  of  Eq.  (83)  is  designed  minimising  the  coupling  effect  ci2.  Fig  14b  shows  the  frequency 
plot  of  the  obtained  coupling  reduction. 
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S12(S)  = 


20600  s2  +52600  s  + 25700 
0.193-10  3  s3  +0.0717  s2  +6.66  s  +  1.78 


(83) 


Phase  (deg) 


£o00=,?ii(s)  (pi*iCs))  1 


L0(s)  =  g22(s)  \P22(S))  1 


Fig.  14  QFT  MIMO  controller  design. 


•  Step  D:  Pre-Filters. 

Open  loop  pre-filters  of  Eq.  (84)  and  Eq.  (85)  are  included  in  order  to  satisfy  time  domain  specifications 
for  reference  tracking. 


/ii(s) 


_ E43 _ 

s2  +  7.5620  s+14.3 


(84) 
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fl2  (s)  “ 


3.026 

s2  +6.355  s  +  3.026 


(85) 


4.4  Experimental  Results 


To  control  the  AdeptOne  SCARA  robot,  a  digital  form  of  the  designed  non-diagonal  QFT  compensator  is 
implemented  in  a  Motorola  68040  microprocessor  (25  MHz),  with  a  8  Mbyte  DRAM  memory  and  a  VME 
bus  card,  and  using  a  sampling  time  of  10  ms.  The  two  plant  outputs,  angles  5,  and  82,  are  measured  by 
encoders,  and  the  two  plant  inputs,  signals  Ui  and  u2,  are  applied  to  a  power  amplifier  that  commands  two 
direct  drives  that  move  the  arms  (see  Fig.  15). 


Fig.  15  Block  diagram  of  the  control  system. 


The  obtained  results  with  the  non-diagonal  QFT  controller,  when  the  reference  ri  for  angle  Si  is 
commanded  from  0  up  to  45  degrees  and  while  the  reference  r2  for  angle  S2  is  kept  constant  (zero  degrees), 
are  shown  in  Fig.  16.  The  same  experiment  with  only  the  pure  diagonal  controller  is  shown  in  Fig.  17. 


Similarly,  the  results  obtained  with  the  non-diagonal  QFT  controller,  when  the  reference  r2  for  angle  S2  is 
commanded  from  0  up  to  45  degrees  and  while  the  reference  ri  for  angle  Si  is  kept  constant  (0  degrees), 
are  shown  in  Fig.  18.  The  same  experiment  with  only  the  pure  diagonal  controller  is  shown  in  Fig.  19. 


Both  real  experimental  results  show  how  the  non-diagonal  controller  (see  Fig.  16  and  18)  reach  a 
significative  reduction  of  the  coupling  effect  with  respect  to  the  performance  of  the  pure  diagonal 
controller  (see  Fig.  17  and  19). 
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Fig.  16  Step  input  at  ri  with  a  fully  populated  (non-diagonal)  matrix  controller. 


Fig.  17  Step  input  at  n  with  a  pure  diagonal  matrix  controller. 
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Fig.  18  Step  input  at  r2  with  a  fully  populated  (non-diagonal)  matrix  controller. 


Fig.  19  Step  input  at  r2  with  a  pure  diagonal  matrix  controller. 
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5.0  QFT  CONTROL  OF  A  WASTEWATER  TREATMENT  PLANT  [1 14] 


One  of  the  main  objectives  of  a  Waste  Water  Treatment  Plant  (WWTP)  is  to  protect  the  water 
environment  from  negative  effects  produced  by  residual  pernicious  substances.  Figures  20  and  21  show 
the  new  activated  sludge  WWTP  of  Crispijana,  Spain,  which  is  able  to  regulate  both  the  ammonia  and 
nitrate  concentration  in  the  effluent,  dealing  with  water  influent  of  about  5000  m3/hour. 

The  control  strategies  designed  to  regulate  that  WWTP  were  based  on  a  hierarchical  structure  where  a 
high-level  or  supervisor  selects  the  set-point  of  the  low-level  or  conventional  controllers.  The  design  of  the 
controllers  was  carried  out  using  the  Quantitative  Feedback  Theory  (QFT). 

Nitrate  control  aims  at  the  optimal  use  of  the  de-nitrification  potential  at  any  moment.  For  this  purpose,  the 
control  algorithm  continuously  adapts  an  internal  recycle  flow  in  order  to  maintain  a  desired  nitrate  set- 
point  in  the  anoxic  zone  (second  loop).  Ammonia  control  aims  at  maintaining  the  required  average 
concentration  of  ammonia  in  the  effluent  by  manipulating  the  Dissolved  Oxygen  set-point  that  commands 
several  air  flow  turbines  (first  loop).  Mobile  average  values  of  some  variables  were  also  introduced  in 
order  to  eliminate  the  perturbations  associated  with  the  daily  24-hours  profiles. 

The  controllers  were  designed  and  verified  using  long-time  dynamic  simulations  based  on  a  multivariable 
and  nonlinear  mathematical  model  (IWA  n°  1)  previously  calibrated  with  real  data  measured  in  the  full- 
scale  WWTP  during  12  months.  The  results  obtained  in  the  regulation  of  the  pilot  plant  show  a  tighter 
control  of  the  effluent  nitrogen  compounds  and  a  significant  reduction  -energy  saving-  of  the  dissolved 
oxygen  demand,  rejecting  the  plant  disturbances  and  insuring  robust  stability. 


Fig.  20  Wastewater  Treatment  Plant  of  Crispijana,  Spain.  (Courtesy  of  AMVISA). 


The  main  control  objective  of  the  first  loop  is  to  guarantee  the  standard  requirements  of  ammonia 
concentration  in  the  plant  effluent,  fixed  on  a  daily  average  SNH  lower  than  2  mg/1.  Accordingly,  the  set 
point  is  fixed  to  1.3  mg/1.  The  maximum  allowed  value  (saturation  limit)  of  the  DO  variable  is  2  mg/1.  In 
Fig.  22  the  daily  average  of  the  effluent  ammonia  concentration  SNH(t )  is  shown  as  a  dashed  line  and  the 
control  input  DO(t)  is  shown  as  a  solid  line,  both  in  mg/1,  which  are  obtained  with  the  Gn(z)  controller, 
under  a  typical  influent  ammonia  load  and  the  temperature  conditions. 
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Fig.  21  Wastewater  plant  diagram  (D-N  configuration). 

Simultaneously,  the  main  control  objective  of  the  second  loop  is  to  guarantee  the  standard  requirements  of 
nitrates  concentration  in  the  plant  effluent,  fixed  on  a  daily  average  SNO  lower  than  2  mg/1.  For  this 
reason,  the  controller  tries  to  strengthen  the  denitrification  process  in  the  D  tank  to  minimize  the  nitrates 
concentration  in  the  plant  effluent.  Thus,  in  the  present  experiment  the  set-point  of  the  nitrates 
concentration  in  the  D  tank  is  fixed  to  0.5  mg/1.  The  maximum  allowed  value  of  the  IR  variable  is  200% 
of  the  design  value  of  the  influent  flow  rate.  In  Fig.  23  the  daily  average  of  the  nitrates  concentration 
SNO  (t)  is  shown  in  the  D  tank  as  a  dashed  line  [mg/1]  and  the  control  input  IR(t)  is  shown  as  a  solid  line 
[per  unit  of  the  influent  flow  rate]  with  the  G^/z)  controller,  under  a  typical  influent  ammonia  load  and 
the  temperature  conditions. 

Figures  22  and  23  show  how  the  system  is  able  to  regulate  within  the  required  specifications  for  both 
loops,  except  in  the  intervals  when  external  disturbances  saturate  the  control  inputs: 

•  First  Loop.  DO  saturated:  [  7~>12,  13->14,  16->17,  34-^35  ]  days. 

•  Second  Loop.  IR  saturated:  [  19->24,  29-^34  ]  days. 


Performance  of  the  First  Loop.  DO  Level  (solid).  Effluent  Ammonia  SNH  (dashed) 


Fig.  22  First  loop  performance  with  Gu(z). 
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Performance  of  the  Second  Loop.  Internal  Recycle  (solid),  Nitrates  SNO  (dashed) 


Fig.  23  Second  loop  performance  with  G22(z). 


In  these  cases,  the  system  runs  in  the  best  possible  conditions,  obtaining  a  good  performance  when  the 
saturation  disappears.  The  control  system  achieves  satisfactory  performance  of  the  effluent  ammonia  and 
nitrates  concentrations  over  the  whole  range  of  operational  conditions.  It  also  obtains  a  notable  reduction 
of  the  running  costs,  minimizing  the  oxygen  supplied  by  the  aeration  system. 


6.0  QFT  CONTROL  OF  A  LARGE  WIND  TURBINE  [124, 125] 

Large  Wind  Turbines  (WT)  present  a  very  complex  multi-objective  control  problem  that  combine  critical 
reliability  issues  with  non-linear  optimization  matters.  Advanced  QFT  robust  control  strategies  have  been 
thoroughly  applied  in  the  design,  development  and  control  of  the  new  real  Wind  Turbines  of  1500  and 
1650  kW,  made  by  M.Torres  company  (Fig.  24).  The  WTs  are  a  variable  speed,  pitch  controlled,  multi¬ 
pole  synchronous  generator  with  two  controlled  IGBT’s  electrical  power  converters  connected  to  the 
stator.  The  main  dimensions  are  about  72  m  of  rotor  diameter  (blades  of  36  m)  and  65  m  of  tower. 


Fig.  24  Multipole  Variable  Speed  Wind  Turbine,  1650  kW.  (Courtesy  of  M.Torres) 
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Fig.  25  Rotor  speed  control  system  block  diagram 


Wind  Speed  -m/s- 


500  505  510  515  520  525  530 


Rotor  Speed  -rpm- 


Time  (sec) 


Time  (sec) 


Control  Pitch  Angle:  ref(k),1(r), 2(b), 3(g)  (deg  blade) 


Control  Pitch  Speed:  1(r),2(b),3(g)  (deg/sec  blade) 


Fig.  26.  QFT  Control  under  very  high  wind  speed  conditions 

The  principal  targets  of  the  more  than  20  loops  of  the  control  system  cover  aspects  such  as  the 
improvement  of  the  maximum  power  efficiency  for  every  wind  speed,  the  attenuation  of  the  transient 
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mechanical  loads  and  fatigue  stresses,  the  reduction  of  the  electrical  harmonics  and  flicker,  and  the 
robustness  against  parameters  variation  with  a  redundant  fault  tolerance  system.  In  addition  some  critical 
problems  arise  in  the  design  of  the  WT  control  system,  such  as  the  difficulty  to  work  safely  with  random 
and  extreme  gusts,  the  complexity  introduced  by  the  strongly  nonlinear,  multivariable  and  time  variable 
mathematical  model  and  the  impossibility  to  have  a  direct  measurement  of  the  wind  speed  experienced  by 
the  turbine,  because  of  the  high  uncertainty  in  the  anemometer  measurement  and  the  strong  influence  of 
the  blades  movement.  These  set  of  motivations  obliged  the  control  engineer  to  get  involved  in  the  design 
of  every  dynamic  element  of  the  wind  turbine  from  the  very  beginning  of  the  project,  and  to  combine 
advanced  control  strategies  such  as  QFT  robust  control  techniques,  adaptive  schemes,  multivariable 
methodology  and  predictive  elements.  The  actual  tests  that  were  carried  out  in  several  Wind  Turbines,  for 
more  than  three  years,  with  the  proposed  QFT  control  methodologies  showed  very  good  behavior  of  the 
WT,  either  in  low,  medium  or  high  winds  and  even  with  the  extreme  30  m/s  case.  Figure  25  shows  a 
simplified  block  diagram  of  the  rotor  speed,  pitch  angle  and  pitch  speed  controllers.  Figure  26  shows  some 
experimental  results  of  the  TWT1650  with  the  QFT  controllers  under  very  high  wind  speed  conditions. 


7.0  NON-DIAGONAL  QFT  CONTROLLER  FOR  A  3X3  INDUSTRIAL  FURNACE  [127] 

This  section  addresses  the  temperature  control  of  a  3-input  (power  supplies)  3-output  (temperature 
sensors)  industrial  furnace  used  to  manufacture  large  composite  pieces  (see  Fig.  27).  Due  to  the 
multivariable  condition  of  the  process,  the  strong  interaction  between  the  three  control  loops  and  the 
presence  of  model  uncertainties,  a  sequential  design  methodology  based  on  Quantitative  Feedback  Theory 
(QFT)  is  proposed  to  design  the  controllers.  The  methodology  derives  a  full  matrix  compensator  that 
improves  reliability,  stability  and  control.  It  not  only  copes  with  furnace  model  uncertainties  but  also 
enhances  the  reference  tracking  and  the  homogeneousness  of  the  composite  piece  temperature  while 
minimizing  the  coupling  effects  among  the  furnace  zones  and  the  operating  costs  (see  Fig.  28). 


Fig.  27.  Industrial  furnace  and  piece  to  be  manufactured  inside.  (Siflexa,  M.Torres-Spain) 
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Response  T2  of  the  3x3  MIMO  furnace 


Response  T3  of  the  3x3  MIMO  furnace 


Fig.  28  Response  of  the  3x3  MIMO  industrial  furnace,  following  a  reference  cure  cycle  and  rejecting  a 
disturbance  at  plant  output  in  the  first  channel  at  t  =  41400  sec. 

(a),  (b)  Ti  and  OPi.  (c),  (d)  T2  and  OP2.  (e),  (f)  T3  and  OP3. 


8.0  CONCLUSIONS 

Since  the  very  first  ideas  suggested  by  Horowitz  in  1959  until  now,  the  Quantitative  Feedback  Theory 
(QFT)  has  been  successfully  applied  to  many  control  systems:  linear  and  non-linear,  stable  and  unstable, 
SISO  and  MIMO,  minimum  and  non-minimum  phase,  with  time-delay,  with  lumped  and  distributed 
parameters,  multi-loop,  etc. 

The  method  searches  for  the  controller  that  guarantees  the  achievement  of  the  required  performance 
specifications  for  every  plant  within  the  existing  model  uncertainty.  QFT  highlights  the  trade-off 
(i quantification )  among  the  simplicity  of  the  controller  structure,  the  minimization  of  the  ‘cost  of 
feedback’,  the  quantified  model  uncertainty  and  the  achievement  of  the  desired  performance  specifications 
at  every  frequency  of  interest. 

The  first  part  of  the  paper  summarized  the  main  concepts  of  the  Quantitative  Feedback  Theory  (QFT)  and 
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presented  a  wide  set  of  references  related  to  the  principal  areas  of  research.  The  second  part  of  the  paper 
introduced  a  method  to  design  non-diagonal  QFT  controllers  for  MIMO  systems.  Finally  the  paper  ends 
presenting  some  real-world  applications  of  the  technique,  carried  out  by  the  author:  an  industrial  SCARA 
robot  manipulator,  a  wastewater  treatment  plant  of  5000  m3/hour,  a  variable  speed  pitch  controlled 
multipolar  wind  turbine  of  1.65  MW  and  an  industrial  furnace  of  40  metres  and  1  MW. 
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